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Abstract 



ri The uniform asymptotic approximation of Green's kernel for the trans- 

O I mission problem of antiplane shear is obtained for domains with small 

inclusions. The remainder estimates are provided. Numerical simulations 
are presented to illustrate the effectiveness of the approach. 



1 Introduction 



Our goal is to obtain a uniform asymptotic approximation of Green's function 
for a transmission problem of antiplane shear in a domain with small inclusions. 
Exact solutions to singularly perturbed problems corresponding to bodies 
^^ with defects are often unavailable. For complicated geometries involved in 

^^ problems of this kind, i.e. domains with multiple small perforations, numer- 

lyv ical algorithms may become incapable of reaching the required accuracy. Also, 

(T^ when the right-hand sides of such boundary value problems are singular, numer- 

^^ ical procedures can suffer from the same deficiencies. In this case, asymptotic 

^~~^ solutions to these problems are desirable. 

J> The approximation of Green's kernels for regularly perturbed problems, 

for the Laplace operator and the biharmonic operator, was first studied by 

Hadamard in [3] . More recently, the question of uniformity for the approxima- 

C^ tions of Green's kernels for boundary value problems in domains with singularly 

and regularly perturbed boundaries, was addressed in [5]. The uniform ap- 
proximations in |,5] were derived using the method of compound asymptotic 
expansions. 

The paper 7J contains the rigorous proofs and remainder estimates for uni- 
form approximations of Green's kernels, given in 9J, for —A in an rt-dimensional 
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domain (n > 2) , with a single small rigid inclusion. Uniform asymptotic formu- 
lae of Green's functions for mixed problems of antiplane shear in domains with 
a small hole or a crack, are discussed in detail in [8] . In [10] , an extension of the 
theory developed in [3 is made to the case of Green's tensors of vector elasticity, 
for an elastic body with a small inclusion. This was followed by |TT] , where uni- 
form asymptotics of Green's kernels for planar and three-dimensional elasticity 
in bodies with multiple rigid inclusions are given. The paper |llj also includes 
analysis of Green's kernels and numerical simulations for antiplane shear and 
plane strain. 

In the present paper, the new feature of the problem tackled is that on the 
small inclusions we prescribe transmission conditions (the continuity of tractions 
and displacements). The inclusions are assumed to be occupied by materials 
which are different from that of the ambient medium. Compared to previous 
expositions into the uniform approximation of Green's kernels in [3[8j[T0j, where 
the kernels are approximated in the bodies containing small holes we also must 
approximate the Green's kernel inside the inclusions. The analysis also brings 
additional boundary layers when the point force is placed inside the inclusion. 

Below, we illustrate one of the main results in this article, for the case when 
the domain has a single inclusion. Let oj^ be a small planar inclusion, occupied 
by a material of shear modulus fj,j , containing the origin O and which is perfectly 
bonded to the rest of the matrix O^ C K^ whose shear modulus is /iQ. Here, e 
is a small positive parameter characterising the normalized size of the inclusion. 
Consider the antiplane shear Green's function N^ for the transmission problem 
inside the perturbed domain flgUuis. We also use N'-^\ as the Neumann function 
in the unperturbed domain (without the inclusion) and i?*^^^ as its regular part. 
We denote by M the Green's function for the transmission problem, in the 
unbounded domain with the scaled inclusion containing the origin. Also let us 
define the vector function V = (I?i,I?2)^, where the components Vj, j — 1,2 
are the dipole fields for the scaled inclusion in the unbounded domain. 

As one of the results we state 

Theorem The approximation of Green's function for the transmission prob- 
lem of antiplane shear in flg^U cu^ C M^, is given by 



iV^x, y) = iV(") (x, y) + AA(£-ix, £- V) + {27^^lor' log(e-i |x - y |) 

(1.1) 
+eV{e-^x) ■ Vxi?(^^'(0, y) + eV{e-'y) ■ VyR<-^\^, O) + ©(e^) 

uniformly for x. y G 51^ U Wg. 

The structure of the article is as follows. In Section[2] we introduce the main 
notations that will be adopted throughout the text and define Green's function 
for the domain containing several inclusions. Section [3] contains the description 
of model fields used to construct the uniform asymptotic approximations of the 
Green's function. In Section |4J we state and prove an estimate for solutions 
to model transmission problems in an unbounded domain with an inclusion. 



Solutions to transmission problems in a domain with several small inclusions 
are studied in Section [5] Results of these sections will aid us in deriving the 
remainder estimate present in the generalization of approximation ( |1.1[ ) to the 
case of multiple inclusions. Asymptotic properties of the boundary layer fields, 
involving the regular part of Green's function for the transmission problem, in 
the unbounded domain with an inclusion are investigated in Sections [6j Then, 
we consider the uniform approximation of Green's function for the transmission 
problem in the domain with small inclusions in Section [Tj and give the for- 
mal algorithm together with the remainder estimates. Finally, in Section [8J we 
demonstrate the effectiveness of our approach and present the numerical simu- 
lations comparing the asymptotic formula in Section [7] with the finite element 
calculations in COMSOL. 

The asymptotic formulae obtained in the sequel are readily applicable to 
numerical simulations. As an example, Figure [T] shows the comparison between 
an asymptotic approximation and COMSOL computation for the modulus of the 
gradient of the regular part of Green's function for the transmission problem 
in a domain with a circular inclusion, for the case when the point force is 
applied outside the inclusion in a planar body. This Figure represents the 
regular part of the displacement field produced by a point force inside a Cast 
Iron disk containing an Aluminum inclusion: Fig. [Tk, shows the computations 
obtained through the formula ( |1.1[ ) when y e O^, while Fig. II|d corresponds to 
the numerical finite element solution produced in COMSOL. The surface plots 
shown are very similar. 
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Figure 1: a) Modulus of the gradient of the regular part of the Green's function, 
computed with the aid of the asymptotic formula ( 1.1 ) when y e 51^, b) A finite 



element computation (produced in COMSOL) for the modulus of the gradient 
of the regular part of the Green's function for the transmission problem in a 
domain with an inclusion. Here, we consider a circular cylinder of radius 30m 
containing a circular inclusion of radius 7m, the shear modulus of the inclusion 
is /x/ = 2.6316 X 10^*'Nm^ (an Aluminum inclusion), where the shear modulus 
for the rest of the matrix is fio = 5.6 x lO^^Nm^^. The position of the unit point 
force is y = (10m, 10m), which is quite close to the boundary of the inclusion. 



2 Main Notations 

Let i7 be a subset of M^, with smooth boundary dil and compact closure Cl. Also 
let oj'^^^ C K^, have smooth boundary 9a;'--''' and compact closure u}^^\ whose 
complement in the infinite plane is Cw^^^ = ]R^\aj(-'\ j — f , . . . ,M. The sets 
uj^'-'', j = 1,...,M, are assumed to contain the origin O, and the maximum 
distance between O and 9a;'-'^ is f. Let ujg , be a subset of fl, with centre 
O*--'-', 1 < J < M. We assume that the minimum distance between O'^-'^ and 
0^''\ k ^ j, k = 1, . . . ,M and the minimum distance between O'^^-' and dVl is 
1. We relate the domain uje' to w^^^ via LOe' — {x : e~-^(x — O'^^) e cj'^^}, 

j — 1,...,M. The perturbed geometry is defined by fig = ri\lj . We , and 
we say that this domain is occupied by a material with shear modulus jiq 
and the domain cje is occupied by a material with shear modulus /i/., where 
(IOt ^J'I > 0, 1 < j < M . In the subsequent sections, along with x and y we will 
also use the scaled variables $,j = £~^(x — O^-'^), rjj = e~^{y — O^^^). 
By XT we mean the characteristic function of the set T, that is 

r 1 , if x e T , 
^ ' \ , otherwise . 

Our goal is to obtain uniform asymptotics for the Green's function N^, of the 
transmission problem in IJ w^ U fl^, which is a solution of 

AioAxiVe(x,y) + 5(x-y)=0, ^ e n, ,y e\J4^ Uil, , (2.1) 

/i,^AxiV,(x,y)+(5(x-y) = 0, x e a>(^) ,i = 1, . . . , M,y e (Jc^Wur!, . (2.2) 

I 

The normal derivative of Ng on the exterior boundary satisfies 

Aio|^(x,y) = -^, xear!,yeUc^Pul], , (2.3) 

where \dil\ is the one-dimensional measure of the set dil; d/dn^^ = n ■ Vx is the 
normal derivative, where n is the unit outward normal. 

Assuming that the small inclusions tOe , j = ^, ■ ■ ■ ,M are perfectly bonded 
to the matrix fl^, we write transmission conditions across the interface doje in 
the form 

ioT j ~ 1, . . . ,M and ye IJ, cje U rj^; the notation duje indicates the exterior 
or interior boundary of the set uJe . 



The symmetry of Ng^, i.e. 

7V,(x,y)=iV,(y,x), 
is guaranteed by the condition 

iV,(x,y)d5x = 0. (2.5) 

The proof of symmetry of A'^, for this problem, is addressed below. 

The symmetry of the Green's function N^ 

Proposition 1 Green's function N^ for the transmission problem in [Ji^i ^e U 
fig satisfies 

7V,(x,y)=iV,(y,x). (2.6) 



dn 



Proof. Let 

r(0) 



TVr^ (x, y) = iV,(x, y) , for x e r!„ y e U;=i ^e U fl 



M ,.,(0 



t(i^), 



,U) 



\M ,,(0 



(2.7) 



m'>i^,y) =. iV,(x,y) , for X e c^^-'',y G UZi^e U n, . 

We apply Green's formula for the functions N^ ' and N^ ^ , j — 1, . . . , M: 
^loj JA^f )(z,x)A,7Vf )(z,y) -iVf )(z,y)A,iVf )(z,x)}dz 

M 

+ Y.^"' ,„ {iVi'')(z,x)A,7Vf')(z,y) -7Vi^')(z,y)AzA^i'''(z,x)}dz 

= ^o/ {^i°)(z,x)^(z,y)-7Vf)(z,y)^(z,x)}d5. 
Jan '- 0"'z oJ^z J 



M 



So;' 



Z=l-'"'^^'' 



M/, 



a^v, 



(/<) 



(9r 



{z,y)- no^, (z,y 



9n, 



-TVi^-: 



lz,yj 






dS'z , 



where the right-hand side is zero as a result of the transmission conditions (2.4) 



the exterior boundary condition (|2.3| and the normalization (|2.5|). Thus 
= /io 



, x) A^iVi^) (z, y) - nP (z, y) A^iVi^) (z, x)} dz 
+ Em/, /^,^ {A^i^')(z,x)A,iVi^')(z,y) -iVf')(z,y)A,iVi^')(z,x)}{2^) 



The next step involves using the governing equations (2.1) and (2.2 1 along with 

the a 

gives 



the above definitions of Tvi''' and N^^'\ j = I, . . . , M. When x,y e n^, (|2.8| 



7Vf'(x,y) = A^f)(y,x) 



(i) 
whereas if x G fl^, y G idf 



iVi^) (x, y) = TVi^^ ) (y , x) , j = 1 , . . . , M 
^e , ye r^e, we deduce 
A^i^^)(x,y)=iVf)(y,x), j = l,...,M 



Similarly, for x G oje , y e fi^, we deduce 



and if x,y G uje 



iVi^^)(x,y) = N^^^y,^) , J = 1, . . . ,M . 
Finally, with x G cj^ , y G w^ , fc ^ j, we have 

iVf^)(x,y)=7Vi^'')(y,x), 1 < j, fc < M, fc^j 



The above relations together with (2.7) lead to (2.6). D 



3 Special solutions in model domains 

The asymptotic algorithm uses special fields defined in model domains including 
the unperturbed set and the exterior of a scaled inclusion. 

1. The regular part i?'^^-' of the Neumann function N^^'^'> in Q is defined as a 
solution of 

MoAxi?('')(x,y)=0, xGf!, 

Qp(r2) a 1 

.o^;^(x,y) = -^((2.)-log|x-y|) + ^, xGa«, (3.1) 

where y G il. 

To guarantee the symmetry of R^ ' we impose the orthogonality condition 

/ R('^){^,y)dS^ = --^ f log|x-y|d5x. 
Jan ^^Mo Jan 

The function N^^^ is related to R^"^ by 

7V(")(x,y) = -(27r^o)-ilog|x-y|-i?(")(x,y). (3.2) 



2. The next model field is Green's function for the transmission problem in the 
domain CHi^^'' U uj'-^\ j = 1, . . . , M. This function is denoted by A/"'^-* and, 
for T] G Cu)'^^^ U uj'-^\ is subject to 



^J■o 



A^Af^^^ (I, v) + S{^-v) = 0, ^ e Cw(^') , 



where the transmission conditions across the interface of the inclusion are 
given as 

on^ on^ ^g^^ 

and at infinity we wiU prescribe the condition 

A^^^'Hl,r,)--(2^Mo)-i log III +c(r,) + 0(111-1), as ||| -> oo . (3.4) 
Symmetry of Af-^K We choose 

civ) = , (3.5) 

and then it can be shown that 

AA(^)(t77)=AA(^")(r,,0, 
i.e. A/'*--'^ is symmetric. 

The proof is analogous to the one of Proposition [l] 
We set 

A/'(^'°) (I, r?) = N^'^ (I, J7) for r? e Ct:^*^) , | G Ccj^^^ U cj^^') , 
AA(J.^) (I, ,7) = A/-0') (I, rj) for 77 e w^^^, | e Cw*^) U J'^i , 

and 

c(0)(j7) = c(j7) for J7eCa'(^') , 

c (■'?) = c{ri) otherwise . 

Assume r/ G Cw^-'^ and let -B_r — {$, : ||| < i?} be a disk, with suffi- 
ciently large radius R, so that q g Bji\il}^^\ By applying Green's formula to 
7V(^'(z,|) and7V(j')(z,J7) in Bi^\Cu^^^ and w'^'' we obtain 

MO / f AAt^^z, I) A,AA(^"'0) (z, r?) - M^^-^^ (z, ,7) A.A^^^) (z, |) | dz 

+M/, / |AA(-'')(z,0A,AA(^'°)(z,J7)-AA(^)(z,»7)AzAA(^''(z,0|rfz 
= ,0 I {AA0-)(z,O^(z,r,)-AA0-^-)(z,r7)^(z,o}rf^.. 



The transmission conditions (3.3) for Af^^^ imply that the integral over du}'^^^ 
is zero. When | e Bfi\ui^^', 



Taking the limit as i? — > oo and employing ( |3.4[ ) , we deduce 






1 51og|z|-i ^ ^/ 1 



which is equivalent to 

AAW(|,r,)-AAW(77,0 = (c(°H^) - c^^HO) lini / j^i{27r)-Hog\z\)dS, 

= c(o)(r7)-c(o)(0 = 0. 

Hence Af^^'^i^,!']) is symmetric for ^,rj G Coj'^^K In a similar way, it can be 
shown that 

A/'(^)(^,r7) - A/'(^)(J7,^) = c^'^rj) - c(°)(6 for I e Cw^^'^ 77 G w^^^ , 

AA(^)(^,r7) - AA(^)(J7,^) = c(°)(r7) - c(^)(0 for V e Ct:^^^^^ G w^^^ , 

Af^^\i,V)-Af'-^Hv,0 = c(^'(^) - c^^HO for € e tj(^\r7 G w^^^ , 



and the condition ( |3.5[ ) implies that the above right-hand sides are zero. Thus 
Af^-'' is symmetric. 

Regular part of Af^^\ Let the function A/''--'-' have the form 
where 



AA(^"^o)(|, r,) = -(2^/.o)-' log II - »7l - /^^^'''^(t »?) 
AA(^-^^) (I, r,) = -(2^A^,J-i log II - 77I - /^i^''^ (I, »7) 

and here h^' ' and ft.]^' are the regular parts of A/'*-^'*^-' and J\f^^'^\ j 



(3.6) 







for|,r;GCw(j) , 

fori G w(j),J7G Clj(J') , 


and 

(3.7) 






for|GCw(^),r;Gw(j') , 
for 1, 7j G U}^^'> . 





Then, the above definitions for /i)^' and h^' lead to this representation 
for 7V(J) 

+X.(.) (l)Xc^u) (»7) ( " 7Z^ log II - 17I - /ii^'''°^ (1, 17) 

+XC..O) (OX.U) (»7) { - ^ log |€ - »7l - /^i^-''''^ (t »?) 

+X.U) (£)X.0) (»?) { - ^ log |€ - ryl - /^i^''-'^ (1, 17) 

The symmetry of A/''-'', then implies the conditions 

h%'°'°^ ii: V) = ^S^'°'°' (^, for ^, ^ G Ca;(^> , 

/i^'''°'(^,r7) = /.S^'°'''(r7,0 + l!-\- ^jlogl^- »7l for r, G Coj^^\ie J^'> 

An { fii. fio J 

and h%-'''^ it V) = h%^'''^ in, i) for ^-q e J^'> . 

(3.8) 

3. We also make use of model solutions known as the dipole fields Vj^' , k = 1,2, 
j = 1, . . . , M, which play the role of the boundary layers in the asymptotic 
algorithm. Let pf-?) = (v[^\v!f^)'^ , where 

T^^'Hi) = Xc^<.)(02?(^'^°ni) +X.<..(I)2?(-''^^H0 , 
and the vector functions V'^^'^^ , I?(J'-f) solve the problem 



/i,^.A^I?(^'^)(0 = O, ieuj^ 



The transmission conditions on the boundary of the inclusion o;'-'^ are 



M/, ^^ (0 - MO ^.^ (0 = (a^/, - /^o)n(^) , on aw(^) , (3.9) 



dne dnc 



v'^''^\i)^v(i^'\i) , onaw(^) , 



where in (3.9), n^^' is the unit normal to w'-^'. At infinity the vector function 
Pfj^o) satisfies 



pU.^)(^) = 0(|||-i), as III ^00. (3.10) 



4. We introduce the function (^^^ which is a solution of 

C(^)(|) = (2^Mo)-'log|^|+C£^+0(|^r^), as 1^1 ^oo, 
where Coo is a constant. 

4 An estimate for solutions to transmission prob- 
lems for antiplane shear in unbounded do- 
mains 

The next result plays an important role in the asymptotic algorithm. It will 
allow us to obtain estimates for the boundary layer fields and derive the estimate 
for a solution of the transmission problem in a domain with multiple inclusions. 

Lemma 1 Let U^^> he a solution of the transmission problem 

U^ (l)||ga„(i)+ = U ■' (Dl^ga^u)- , 

U^^\^)^0 as lll^oo, 

where tp^^' G L^{duj^^'), d/dn is the normal derivative on the smooth boundary 
du!^^', outward with respect to uj^^\ and 






(4.1) 



We also assume that 

where Q^' is given as solution of Problem 4 o/ iSecfionpl Then 

sup {(1^1 + l)\U^'\m < const ||<^(^)L^(a.o)) , (4.2) 

where the constant depends on fio, l^i o.'^T'd du^^', for j = 1, . . . , M . 
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Proof. We note that transmission problems are studied in detail in [T] and 
[2] , in the context of boundary integral equations and their solvability. 

Let us first represent the solution £/*■•'■' by two functions [/(■'•'^) and U^^'^\ 
harmonic in the domains Cui^^' and oj^^', respectively. These functions satisfy 

For the function U^^''~'\ the condition 

[/(j,o)(^)^Q g^g lll^oo, (4.3) 

also holds. 

Applying Green's formula to the functions A/"^^-* (see Problem 2, Section^, 
t/(J''='' and [/(J'-f), one obtains 

UU.o) (!) = _/■ AA(^) (C, »7)^^^Hr7) d^r, , (4.4) 



for^e Cw(^), and 



t/0-.^)(0 = - / AAW(l,r,)^(-')(r,)d^^ 



(4.5) 



for ^ e w'-'^. 

First, let |^| > 2, then using the asymptotics for the function A/"^^-* at infinity, 



and the condition (4.1), we deduce 



(1 + |||)|C/(^'°)(0I < const ((1 + |C|)| log llll / V^^'Hrt^dSr, 

+ll<^*^^llLi(ea.o))) 
< const\\ip'-J^\ r(oujuy) . (4.6) 



Also by the Cauchy-Schwarz inequality and (4.4) 



|C/(^'O)(0| < const ||¥.(^')|L,(a.u,) < const||^(^")|L^(a^u,) , for ^ G B^'^ , 

(4.7) 
where B3 = {^ : ||| < 3}. 



Similarly the integral representation (4.5) gives 



|[/(j'-f)(^)| < const IML^^Q^ur) < const\\ip'^^^\\^^^g^uy) , for ^ G w 



(j) 



The combination of (4.6), (4.7) and (4.8), leads to (4.2) 



(4.8) 

n 



As an immediate corollary of Lemma [Tl we have an estimate for the dipole 
fields associated with the scaled inclusion uj^^\ j — 1, . . . ,M: 
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Lemma 2 For the dipole fields V ^ , j — 1,2, i = 1, . . . , M , 
sup {(1^1 + l)|Df(0|}< const , 

holds, where the constant in the right-hand side can depend on jiq, fJ-ii and oj^'^' . 



5 An estimate for the maximum modulus of so- 
lutions to transmission problems for antiplane 
shear in a domain with several small inclu- 
sions 

Here we obtain an estimate for solutions to transmission problems for antiplane 
shear in domains with small inclusions. The next lemma will be used in Section 
[T] incorporating the remainder estimates produced by the approximations of 
Green's function in a perforated domain. 

Lemma 3 Let u be a function in Loo(Ui=i'^e ^ ^e) such that Vu is square 
integrable in a neighbourhood of doji , i = 1, . . . , M . Also, let u be a solution 
of the transmission problem 

fj.oAu{x) = , X e Oe , 

fii^Au{x) = 0, x£uj^J-\i^l,...,M , 

Qu 

MOT^(x) =^(x) , xedn, 



du : 



du . 



- m(0 



w(x)| ^fl (0+ — u(x.)\ ^„ 



V'e\^) 



(5.1) 



where ^ g L^idQ), if I' G L^idtj^'), fori <i<M, 

[ v(x)d^x = o, / 4''>{^)ds^^o 

Jan J9w<'' 



and ipj' (x) — e ^Lp^'^'[e "'"(x— O*-*')), i — 1,...,M. To provide uniqueness we 
also assume 

I u(x)d5x = . 
Jdn 

Then there exists a positive constant A, independent of e and such that 



ll"llL^(Ufii4''ua,) ^^jll^lli=c(S")+S^'^^M"'^^ "i=o(ac.f')f ■ (^-2) 
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Proof. a) The inverse operators to model problems in fl and Cuj^^' , j 
1, . . . , M. Let us introduce the operators 



N -.i^^w and m^''> : ip^^'^ - 
which are the inverse operators of the problems 

^o^'w{x.) — , X e ri , 

dw , 



(5.3) 



and 



/"o^(x) = V'(x) , xedn, 

[ w(x)dS'x = 0, 






(5.4) 









^'' dn ^^'kedujU)- 



- m(j) 



(pU)(|), 



^eaLj(3)+ 



„UJ(|) 



$e9w«)- ' 



(5.5) 



t;(j)(O^0 as lll^oo, 
where ip & Lao{dVl), Lp e ioo(9a;(-''), j = 1, . . . , A/, also 

/ ^(^'^d) dS-g = and / V(x) dS-x = . 
JSwO) Jan 

In scaled coordinates ^, = e^^(x— O'^-'^), j = 1, . . . , M, the operator Die is 
defined by 

(9TW^F)(x) = (^(^■V^^'')(«,) , 

where pi^ (y^) = £-V*^'H£"^(x- O^^'))). 

6j j4n estimate for solutions to the model Neumann problem in fl. Let 
iV^^^ denote the Neumann function (3.2) in fl. 

Then an application of Green's formula to A^'^^'(x,y) and w{y) yields the 
representation for w 



u;(x)= / 7V(")(y,x)^(y)d5y + 



1 



w{y) dSy . 



(5.6) 



The solution w of the Neumann problem in fi is subject to the orthogonality 



(5.6) is zero, so that 



condition in problem (5.4), and hence the last term on the right-hand side of 



an 
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From this we obtain the estimate 



sup \w\ < const sup 1-01 



(5.7) 



c) The case of the homogeneous boundary condition on d^. When the 
right-hand side of the Neumann condition on dfl in (5.1 1 is zero we look for a 
solution of the form 



M 



M 



d 



Ui 



(5.8) 



where g^ is an unknown function defined on duje such that 

/ g<^^\^^)dS^^^O, 

JdijU) 

and g^-'\^j) = egP (x). The function ui is harmonic inside IJi=i "^e ^ fi^, and 
is continuous across the boundaries of the small inclusions. On dil, ui satisfies 

/^o^— (x) =0 . 
on 



Computing the jump in tractions of ui on duJe , we obtain 



dui 



dui 



where g^x) = (g^^H^), ■ ■ ■,g^^''H^)f and 



5^) + (5i™)5s)(x) 



(5(")g,)(x) = (^o-/^/„j[^{ E "^^9^^} 



Ix^duj^ 



d_ 
dn 



l<j<M 

M 



{N(TraoMoEa;^(^i^'^3^^^))}lx, 



j=i 



(5.9) 



Let iJ^™) denote a disk centered at O^™^ and containing tOe ■ Using a 
local estimate for solutions of Laplace's equation, along with Lemma [l] and the 



definition of ge above, we obtain 
d 



dn\ 



E "^'al^^} 



l<j<M 



LooCa^"') 



< const II E ^e^9i 



U) 



l<j<M 



< const 2. £ 



'U'^^ 



e IlL. 



,r.,„«).(5-10) 



l<j<M 
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Then, from the local estimate for harmonic functions, we can also assert the 
inequality 



dn 



{N(Traa/.oEa„(^i^^5i'^))} 



M 



< const \\N[Tran Mo Y. a^I^^^'^^i^'^: 



j=i 



M 



^ const 5] -(^(^-'ff^^); 



j=i 



_9_ 
I (9n' 



L„„(a4'"') 



ioo(B('")) 



ioo(an) 



< const E-^'ll3e''^llL„(a<.«') 



where in moving to the second inequality we used estimate (5.7), and then 
Lemma [T] brings us to the last inequality. 



Then, the preceding estimate and (|5J|), ( |5.10[ ) lead to 

II ^i")-?. I! 



M 



'L^idujl"^') 



(„), < const e^^ II "(J) 



5s 'llL^(a4^)) 



j=i 



Hence, from the smallness of S'e , we can write 



(m) 



ge = (I + S,)-Ve 



.(^^)i^\\T 



where 5s(x) = (gi^^(x), . . . ,5f^^(x))^, ^.(x) = (^i'^(x), . . . , ^^^^(x)^, S, is 



(1) 



JM) 



a matrix whose rows are Se , ■ ■ ■ , Se , and 



i(j) I 



ffr'^llr /fl (j)^ < const max a? 

^e "L^(du:y') - l<k<M" 



(fe)|| 

'^ Ilio 



(5.11) 



From (5.8), together with (5.7) and Lemmalll we obtain 



"il 



i=o(Uf£i4''uo,) 



M 

< const Y.^\\9i'^hud^i^') ■ 

3 = 1 



Now, this and (5.11) give 



\ui\\. ., ,M (0, ,r, ^ < const £ max llcoJ.'^'lL ,„ (»=), 



Sj) 



(5.12) 



d) The case of continuous tractions on duie , j — 1,...,M. In this 
situation we look for the solution U2 in the form 

U2^'NiP + v. (5.13) 



Then, z; is a solution of (5.1 ) with the boundary conditions 



dv 
/io ^ (x) = , on dn , 
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and 



i;(x)| 



ceOc 



= w(x)| 



cGa< 



.0)- 



MO 






for 1 < j < M ^ where the right-hand sides of the above traction conditions on 
9we , j = 1, . . . ,M are self balanced, so by part c) of the present proof 



Mr i\\^' ,,(''iio ^ ^ const e max \\^- — (x) 



ioo(94'°') 



< const e max Nf/^L ,„,i.,, 

< const e\\ip\\L^{dn) ■ 



This inequality, (5.7) and (5.13), give 



ll"2llL^(Ufli4"uO,) ^ ^^O^^t ll^lUooOfi) • 



Finally, we obtain (5.2) through the combination of (5.12) and (5.14). 



(5.14) 

D 



6 An asymptotic approximation for the regular 
part of the Green's function J\f^^^ at infinity 



In this section, we shall prove a result which concerns the asymptotics of /i^ 

, '^) a — 1 ]\/r tcr^r^ /Iq rV\ n+ ;„flv,;f,r 

'at 



0,0) 



and n"^' ' , j = 1,. . . , M, (see (3.6|), at infinity. 



Lemma 4 For |^| > 2. Then 
a) 

h^-''\tv) = -2?(^'°'(r7) • V«((27rMo)-Mog|iri) + 0(111-2(1,71 + l)-i) 



for T] e Ci:^^-'), 
b) 



'■JV 



(l,r,) 



-(m7/ - Mo'){(2^)^' log 1^1 - »7 • V«((2^)-i log III)} 
-2?(^"-^) (r,).V«((27rA.o)-' log lir^) + 0(111-2). 



/or ri <E uj 



U) 



Proof. First we study the asymptotic behaviour of the functions h^' ' , 

/ijv' ' : 'ijv' ' introduced in Problem 2 of Section K3l in the neighbourhood of 
infinity. We show that for ||| > 2, 
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a) 

h^i-''-''\i,v) = -P(^'°'(r7) • V«((27r^o)-'log|€r^) + 0(|^r^(|r7| + 1)-^) , 

(6.1) 
for rj e Cw^-'), 
b) 

h^''-''\v,^) = -V^'^'Hv) ■ ^di^nf^o)-' log m-') + Om-') , (6.2) 

for T] G a;'^-'^ 
c) 

h'-A-'^-'Htv) = -(M7/-Mo'){(2^)''logl^|-»7-V«((2^)-Mog|||)} 

-P(^"'^)(r,) . V«((27rA.o)-^ log lir^) + 0^-') , (6.3) 

for rj G o;'^-'^. 

i) Auxiliary functions h^ and T^^'. Let us introduce the function h^ as 

From Problem 2 of Section K5l when ry G Cw^-'-' U w'^^'-', /i];J is a solution of the 
transmission problem 

/^^.^;^(^'»?)|«eac.(.)--/^o^;^(l,r7)|^^g^,,,+ 



27r V A^o M/j / <9"C 

^i^'(^'^)lcea.o)+ = 4^(^'^)l€Ea.o)- ' (6-6) 

/^i^^ (I, r,) = -(M7/-Mo')x.u)(»7)(2^)-Mog 1^1 +0(111-1), as ||| ^ (^ . 

We introduce one more auxiliary vector function T'^-''(^) — {T^ (^)}|^j^, 
defined by 

T^'\i)=i-V^'\i). (6.7) 
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It solves the transmission problem 



' dn. 






(6.8) 



T(^)(0=l + O(|ir^), as lll^cx), 
which is consistent with Problem 3 of Section [3J One can also represent T^^\ 

a) The asymptotics of h^' ' at infinity. For |^| > 2, r) £ Cuj^^\ we 



as 



have from Lemma 2 in [6], that the function h^' ' defined in (6.4), has the 
asymptotic representation 



Aj,o,o) 

''N 



itv) = C(^"^°)(r7) . V«((27r)-i log \^\-') + 4' ' \^^ V) , (6.9) 



where the remainder r]:!T'^'^'{^,r]) satisfies 



N 



S3,0,0) 

' N 



(|,J7)| < const (1 + \r]\)-^\^\-^ for ||| > 2,7? e Clu^J^ 



The vector function C''-'''^'(j7) = {Cj (^)}f=i is evaluated below. 

Evaluation of C^^''-'\q). Let Bji be a disk centered at the origin with a 

(i) (1) 

sufficiently large radius R. We apply Green's formula to /i]^ and T^ ' in the 

domain Bji\lj'-^'> U oj^^-* 

= mo/ {Tp^(|)A^4)(tr7)-4'(er7)A4Tp^(|)}d^ 

+M/, / {T(^''(OAc/ii^^(^,r,)-/^i^)(|,77)A«T,p)(0}d| 






dh 



U) 



dh 



U) 



^^'^^ \ ^^^ in^^^' '^^l^^"'- - ^«in7(^' ^)|^^«' 






dS^ , (6.10) 
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where while combining the integrals over dio^^^~^ and du)^^'^ we have used the 



continuity conditions for the functions h^ and T^"* , respectively (see (6.6) and 



problem (6.8 1). With the use of the transmission conditions for T\ ' and h 



U) 



when rj £ C(jJ^^\ this equality becomes 
= Mo 



N ' 



/».{"'• ''«?!'«••" -'■~'«^''>?l^'«}'*' 



(a*/, - Mo) 



Mo J9w(j) 



d_ 



T^P{i)i^{{27:)-'\og\S,-'n\)dS^ , (6.11) 



for r) e Cuj'^^y The last integral, owing to (6.7), is equal to 

= ^^^^^|6^((27r)-Mog|^-r7|)-I?p)(0^((27r)-Mog|^-r7|)|d5« 

where rj S CH)^^'. Now, using the jump in tractions for the dipole fields, stated 
in Problem 3 of Section |3j we can write 



f^O Jdu 



1 



^J■o Jdu(j) 



(27r)-ilog|^-r,ri 



M/,^^(0-MO^^(I) 



iO)/ 



d 



-(M/, - Mo)I?r^(0^((27r)-i log II - v\-')j dS^ . 

Using the continuity of the displacements for the dipole fields (see Problem 3, 
Sectional, we write the above right-hand side as 



M/, 



(2^)-Mogi^-7?, 



i5^^''(^)_I,p./)(^) " ((2.)-Mog|^-„|-)U5, 



_9_ 



-/ i(27r)-Mog|^-r,r 
Jauj(Ji [ 



im 



U,o) 



d 



dn^ (0-I'?'°'(0^((2^)-Mog|C-r,r^)^dS« 



Then, upon applying Green's formula to P, , V^' in Bb\ui^-^\ uj^^\ respec- 
tively, with the fundamental solution of —A, and using the definition of T'^-'^ 
we obtain the above in the form 



V\'^''\'n)~ I |(2^)-Mog||-r,| 



-idV, 



U,o) 



dnc 



-(0 



-v\''''\i)^ii'^^r'^og\i-'n\-^)]ds^ 
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where ry G Cui^^^ and by (3.101 the integral over dBj^ decays as i? — > oo like 



0{\\ogR\/R). Therefore, we have 
(a^/j - Mo) 



Mo 






_ T>{3,0) 



iv) 



for T] £ Cijj^^\ Equation (6.11), with substitution of the preceding equality, 
leads to 



I>.'*°'(») = w 









dn^ 



(6.13) 
We now aim to determine the vector function C^^'^^q) — {C^ {'n)}^^i in 



(6.9). Taking the limit _R — > oo in (6.13) and using (6.9) (where h^' ' (C^?) 



h^it ri) for I, J7 e Ccj(j')), we have 



-Pp'°)(,7) - lim MoC(^^o)(r/)- / {e.T^(V«((2^)-i log l^r^)) 
-V^((27r)-Mog|^ri)^}d^^ 



= MoCp'°)(r,) 
Thus, we have derived that 



C0^o)(^) = _ 



2?0-,o)(^) 



Mo 



which corresponds to the leading order term of the function h^' ' , stated in 



the current lemma in (6.1 ). 



Hi) The asymptotics of h^i-^-'^i at infinity. We have already shown that for 
III > 2, J7 G Ca'(j') 



Ao,o,o) 

'-N 



it V) = -V^^'^'Kn) ■ V«((27rMo)-^ log lir^) + r^i^"^"\i, ri) . (6.14) 

In order t o de duce the leading order term of h'^J^^^'-'^ at infinity, we recall the 

is symmetric for ^,ri £ Cuj^^' we have the 



relations (3.8). First, since h 



U,0,0) 

N 



A],o,o), 



above asymptotic representation also holds for h^^'^'^ [r],^). Next we allow ry 
to approach the boundary of the inclusion oj^^K For rj £ dui'^^^ we have 



I?(j'^o)(j7) ^V'^3J)(^r,) 



'A],o,o) 

'■N 



iv,0 






(»7,l)- 



Therefore, allowing 77 £ u)'^^\ for ||| > 2, we arrive at 



N 



itv), (6.15) 
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N 



is the remainder term and subject to its smallness the leading 



where r 

order part of (|6.2|) has been formally deduced. 



Remainder estimate. Consider the function 
which by ( |6.14| and ( |6.15| ) is 



,0) 

' N 



itv) 



+X.U) (r/){/ii^'''°'(r/, + 2?(^'^)(r,) • V«((27rMo)-^ log l^r^)} . 



Let 1^1 > 2 and write the problem for rj^"^ with respect to tj as follows 



U) 



IJ-o^ 



Ar,4\tri)^0, v^CCo^'\ 



So) 



Hi^Ar,r'^>i^,ri) = 0, rye 



,{j) 



dr 



U) 



^'' dUr, ('^''')|^69-"' 



Mo 



9r 



9n^ 






ri€du:(3) + 



= ^^'' ^"U ^((27r)-4og|g-.7|)-A((2.)-iiog|g|)},(6.16) 
/iO I- wrif one J 



„(i) 



0)/ 



' N 



^^'''^>\r,eauJ(l)+ ^"^N (^''7)|T,Gaw(3 






(^r,) 



O, 



|J7| — ^ oo 



We have that the right-hand side of (|6.16|) is self balanced and 

d 



(m/, - Mo) 



Mo 



|A((2^)-ilog|^-r,|)_-^((27rrMog||| 



< const \ri\\$\ ■^ < const |^| "^ 



where -q e d(jj''^\ \ri\ < 1 and |^| > 2. Now an application of Lemmalll leads to 



the estimate for r 



|rj^^(l,r7)|< const 1^1 



1^1 



for 1^1 >2, r;e Cw^J) U w^^). 



iuj The asym,ptotics of h^' ' . We once again refer to ( 3 

2tt [fii^ fio 
For 1^1 > 2, r) e a;'-'', this can be rewritten as 



for the relation 



'TV 



h^f;r'^^"itv) = h 



N 



log II-T7I for ^ G Clj(^), r/ G w^^') 



Ai.oj) 

''N 



itv) 



AjJ^o), 



h\'.--'(n,i)- 



-I- 



^}{log|l|-r?-V«(log|||)}+(^^ 



By combining this with (|6.2| we obtain (6.3) 



D 



The proof of Lemma I4] is then completed by applying (3.7), (6.1) and (6.3). 
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7 Uniform approximation of N^ in a domain with 
multiple inclusions 

The aim of the current section is to present the uniform asymptotic approxima- 
tion for N^ in a domain with several inclusions. 

Theorem 1 The approximation of Green's function for the transmission prob- 
lem of antiplane shear in IJ; We U fig C ffi^, is given by 



iV,(x,y) = iV(")(x,y) + E-^*''(^.''?.) + Mi27Vf,o)-' log(e-> - yl) 

(7.1) 

M 

+e^{D'^'(^^.)-Vxi?<">(0'^\y)+©'^'(^,)-Vyi?<">(x,0''')}+0(e') 



uniformly for x, y G IJ^ w^ U f^e . 

Proof Formal asymptotic algorithm. First, we give a plausible argument 
concerning the representation of A^^. We propose the function Ng to be given 
in the form 

iVe(x,y)^-;lf^^^^+f MlM')iog|x-y|-i?.(x,y), (7.2) 

where for y e U;=i '^e ^ ^ej ^e is a solution of 

AioAxi?e(x,y) = 0, xe^!, , (7.3) 

M/,Axi?e(x,y) = 0, xeco«,i = l,...,M, (7.4) 

dRe, , MO f xnAy) , ^ M'^^yA g(log|x-y|) , 1 ^_ 

(7.5) 
and for i = 1, . . . , M , R^ satisfies the transmission conditions 

l^h-l^o ( xjxAy) I y X4"(y) \ 9(log|x-y|) 

^^(^'y)lxe9^w+ == ^^(^'y)lxe9c.«- • (7-7) 

In the above problem we have that i?g is subject to the orthogonality condition 

/ i?.(x,y)d5.^-^f^;^^+^MlM') / log|x-y|ci5.. (7.8) 
Jan 2tt y no ^ M/i / Jan 
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First order approximation for R^ . If we allow y to be located inside one of 
the M inclusions, then we assume y £ ijJe , where m is fixed, 1 <m < M. We 
first rewrite the boundary conditions for R^. 



The boundary condition for R^ when x G dVt, y G O^ U a;^ , is equivalent 



to 






+x„(-)(y)<^- 



27r duy, 



\ log |x - y| H log |x - y| 



\dVL\ 



Using scaled variables we can also rewrite the transmission conditions (7.6) on 



+X4-)(y){ 



dR, I 

= -xa,(y) 

M/, - Mo 9 
27r duy. 



9^"- -Mo^(x,y)|^gg^(,)+ 

M/, -Mo^(logl^g-^gl) 
27r/io 9nx 

) log |x - y| + — log \i - J7 I 

M/„. Mo/ t^o 



}. 



where y G ilg U w^ . Finally, the transmission condition (|7.6|) on the m 
inclusion becomes 



^th 



dRe, ,1 dRe 

a;^(^'y)lx,a4'"'--^09;i; 

t^im-t^o XnAy) , ^^^"•'(y)^ 9(log|lm-»7ml) 



M/™^;^(^'y)|xea4'"'- " ^o^Jf ^'''^^l-e9c.(")+ 



27r 



Mo 



M/„ 



(9?- 



,(™) 



for y G r^s U oje . We therefore have the representation 



i?e(x,y)=Xa.(y){i?^'^Hx,y) + ^4^'°)(C„r,,)} 

fe=i 

+ E 4'''^(^.,r,,)}+^s(x,y). 

k^7n 
l<k<M 



(7.9) 



r/ie construction of boundary layer terms. The function Z^, given in (7.9), 
is harmonic for x G lJi=i i^e U fie, y G ilg U We and is continuous across the 
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boundaries of the small inclusions. The normal derivative of Z^ on the exterior 
boundary is 



dZ, 



\^o -^ (x, y) = -MoXo, (y) X! 



M p,Ak,0) 



dh 



N 



-Aiox^(™)(y) 



fc=i 



dn^ 



-i^k^Vk) 



a r 1 / 1 



9?- 



1 



271" VM/™ Mo 



log |x - y| 



+ h''^'^\^m,VrrJ+ Yl ^n''^'' i^k, Vk 



k^m 
l<k<M 

for X G dQ, and the traction transmission conditions for Z^ on the small inclu- 
sions are 

dZ I dZ I 

= -(/i/„-Mo)(xn.(y) + X4-.,(y))^{i?(")(x,y)+ ^ 4''''^(€fc,^/c)} 



and 



k^m 
l<k<M 



dZ, 



dZ, 



^^h^^i^^y)L^au,i-^- -A*oi^— (x,y)L ^„ (,) 



'dr 



dr. 



IxGaLjg 



-K-A*o){x...(y)^[i?("nx,y)+ E h''N''\ik.'nu 

k^q 
l<k<M 

^ ^^ ^ ^ ^ \og|x-y|+i?(^)(x,y) 



+x,„(™) (y) 



9nx 



27r \fii^ A*o 



-^i^''"^^(€m,^m)+ 51 ^7v'°''(^fc,^fc) } 



k^q 
l<k<M 



where l<q<M,qy^m and y e il^ U Wg . By expanding the first order 
derivatives of the logarithmic term near x <™) (y), about y = O^™' up to 0{e^), 

in the exterior boundary condition and the traction condition on doje , q y^ m, 
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we obtain 



Mo 1^ — (x, y) = -fJ-oXn, (y) 2^ 



-MoX,„(-) 



(y){ 



dn^ 



k=l 

1/1 1 



dr 



i^kjVk 



1ok|x-0(")| 



.2tt V/x/^ Ho, 
(y-0(")).(x-OM) ^ ( ) (,,o. ^l^n^20 



l<k<M 



for X G cTi, and 



dZ, 



dZ, 






l<fe<A/ 



+X4"^)(y){ 



dn^ 



1/1 1 



27r VM/,„ Mo 



(log|x-0 



_ oMi _ (y-o<"))-(x-o('")) 

|x-0(™)|2 



-i?(")(x,y) + 4"^^)(|„,^„J+ ^ 4-''''(l..,r7,)]+0(e2)} 



k^q 
l<k<M 



where y G ft^U We " . Next we apply Lemma kl in order to rewrite the above 
boundary condition on dfl for Z^ as 



dZ 
Mo^ — (x,y) = -fJ-oXn. 






d reVC'-O^i'nk) x-OC^) 



fe=i 



Af 



Mo 



2Trno |x-0(fc)|2 



fe=l 



9?- 



27r^o |x-0(™)|2 



-(Ai/i-A*o^)log'^+ Yl 



k^m 
l<k<M 



2^T^lo |x-OW|2 



O 



(£^)} 



(m) 



for X e 3fi, y G r^e U Wj . The same lemma in combination with the Taylor 
expansion about x = O*^™' of the first order derivatives for the function R'^^^ 
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leads to 

dz I dz I 

= -(Xo.(y) + X„(™)(y))(M/„ - Mo)n(") • Vxi?(^)(0("\y) + 0{e) , y G f], U c.(") 

where n'-'"-' is the unit outward normal to the inclusion Wg . Similarly, on the 
g*'' inclusion, 1 < q < M, q y^ m, we derive 

= -(xo.(y)+X4™)(y))K-A^o)n(«^-Vxi?('^)(0(^),y)+0(e), y e ^i, U c.(™) . 
Then, we approximate Zg by 



Z,(x,y)^-(xn.(y)+X^(^,(y)) E e{p('=^°)(77fc)-Vyi?(^^)(x,0 



ik)\ 



l<k<M 



+to.(y) + X4™,(y))2?("HO • Vxi?(^^)(o("),y)} 

(7.10) 
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Combined formula for N^. The substitution of (7.9), (7.101 into (7.2), for 



y e r^e U Wj , yields 



M 



iVe(x, y) = xn, (y) -(27^^ior' log |x - y| - i?(")(x, y) - ^ ^'^^ (^„ r,,) 
I fc=i 

M ^ 



fe=i 



+X„("^)(y) -{27rf^jJ-Hog\^-y\-R(^\^,y)~hf''\^^,v„,) 



l<fc<M 

+£p(™'^)(r7„J-Vyi?(f')(x,0("))+e ^ ^'^'''(■^fc) ' Vxi?^''^(0('=),y) 

1<A;<M 



+2?(^'0)(»7fc) • Vyi?(^^)(x,OW) +r(™)(x,y) 



Then, by the definitions of N^^^\ J\f(-^^o) ,^^^ j^(in,i) ^ ^^A^ jg equivalent to 



(7.11) 



M 



I fe=i 

M "j 



k=l 
M ( 



l<fe<M 



M 



-M(27r/io)-Mog(e-i|x-y|) + £E^^'^(^'c)-Vxi?("Ho('=),y) 



fc=i 



£l?(^'^)(r,^-)-Vyi?(^^)(x,0(^))+e E ^^''°Hr7fe)-Vyi?(^^)(x,0(^))|+r,(x,y) 



l<fe<M 
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Remainder estimates for the approximation of N^ . We represent r^ as 

M 

r,(x,y) =Xf2.(yme(x,y) + ^X4p)(y)f)i^'nx,y) , (7.12) 

where SJlg and l)e,j = 1, . . . , M are defined below. In what follows, we estimate 
9Jle, and ije , j — I, . . . , M, in order to estimate r^. 

Remainder estimate for the function 9Jl£(x, y). First, let y G il^. According 



to (7.11), the fimction DJl^ is a solution of 

/ioAx»le(x,y) = 0, xerJe, 

^/, Axmte(x, y) = 0, x e 4'^ * = 1, . . . , M , 
and on the exterior contour is subject to 
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(7.13) 



with transmission conditions on duje of the form 

l<j<M 
j=l ^ ^ X J 

-(m/. - /.o)e5]l?(^'°nr7,) • Vy^^(x, 0(^-)) , 
j=i ^ 

9Jte(x,y)| ^„ (i)- =9Jtj(x,y)| ^„ (i)+ , 
for 1 < i < M, where y G i^e and 



(7.14) 



/ »t,(x,y)rfS'x = 0. (7.15) 

Jan 



I an 
Before estimating the boundary conditions we observe that 



(7.16) 



/ Mo^ — (x,y)d6x = 0, 
Jan on^ 

L^'^ (^^•^("'y^lxea^''- - t^o^i^,y)l,,^.n) dS^ = ^7,.17) 
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forz = 1,...,M. 



Estimation of the right-hand side o/( |7.13 1 ondVt. Since x G Sfl, |x— O'^'j > 
1 for j = 1, . . . , M , and by Lemma [2J the estimate 



M 



^Mo|E^?^(^.)-^xi?(^''(0(^),y) 



j=i 



c^n-x 



< const e , 



holds for X e 951, y G fig. Using Lemma Ilk), one has for x e dfl, 



l^o 



fJ-o 









c)nxV|x-0(j)|2. 
< const e^lx- o(^')|-2(|y _ o(j')| +e)-i < const e^dy - O*^)] + e)~^ 



where we have also made use of the boundary condition (3.1) for R^^^K The 
previous two inequalities then give 



MO 



dm. 



duy, 



(x,y) 



< const e^ X e 917, y e fig 



(7.18) 



Estimate for the right-hand side of (7.14) on due , i ~ 1,...,M. The 



,(») ,■ - 



regular part 7?^^^ of N'^^^ is smooth for x, y G fJ, and we can expand the first 
order derivatives of this function about the centre of the small inclusion loI to 
obtain 



(M/i - i^o)'^j^ (x, y) - e| 



A*/, 



^Mo- 



ao('-o) 



(^. 



,(0 + 



dn. 

}-Vxi?(")(0«,y 



^^i^lxeSt. 



(0- 



9nx ^^''Ixe9t>;j 



(/ii. -/io)n^''-(Vxi?^"'(x,y)-Vxi?^"^(0^'',y))| 
< const £ X G dujy' , y G r^e . 



(7.19) 



With the use of Lemma [2] we have 



E 



i<i<M 



l^h- 



(^j)Lpfl,.,(')- -MO 



9D(^-0) 



(^. 



9nx ^"^^ixeeu.r'- ''^ ^rix ^''^■^ixea4'>+ 
< const £^, for x G duj^^' ,y G f^e , 



Vxii^"'(0^^\y) 



(7.20) 
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and the same lemma also gives 

M 






M 

< const ^e2(|y-0(^')|+e)"\ ior x e duj^'\y e n, . (7.21) 



Due to Lemma |4k), 



(j,o) 






l<j<M 



l<j<M 

Therefore, this estimate with ( |7.19[ ), ( |7.20p and ( |7.21[ ) lead to 



M/, 



9971^ 



-(^'y)lxea4"--Ai0 7^(^'y)lxea4"+ 



< const £, X e (9w^*'',i = 1, . . . ,M, y e fie 



Then, by Lemmalsland the preceding estimate with (7.15|-(7.18l, we obtain 

M 

|SH,(x,y)| <const£2, x e [juj^ U n,,y e ft, . (7.22) 

1=1 

Remainder estima te fo r i)s , j — ^, ■ ■ ■ ,M. Let y e cj^ , where j is fixed, 
I Y (j) 

1 < J < M. Then by (7.1 1, the remainder term f)e solves 

^oAx#(x,y) = 0, xef!, , 

Ai,,A,0(^)(x,y) = 0, X e 4'),z = 1, . . . ,M , 
and satisfies on the exterior boundary 






E ^«....)--^(y-v,«'-(o-., 



\<k<M 



9nx 



a _ „,o., „... ^ r^De^'^^^^^^.^^^coj^ow^y) 



+el''^-^'(r/,) • ^V,i?("'(x,0(^)) +e ^ 



k^3 
l<k<M 



9nx 



+0"='°)(„J.^V,i?(")(x,0(^))|[, xG9f7 



(7.23) 
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Also, on the boundary of the i*'' mclusion, [)e is subject to 






f;^(x,y)-.i^'-'(..)-^v.i.(")(x,oO)) 



9 



e E I'"''°'(»7j-^V,7?(")(x,0(^))+ E (1-M^^^(^.,r,, 



dh 



(k,0) 



l<fc<A/ 



l<fe<Af 



-fe(27rMo)-'^(log(e-^|x-yl)))-(l-5,,)^^(^^,r,^.) 






an, vs,;ixe94"- ^" 9^^ ^''^■^Ixes-i* 






„(•)+ 



Vi?^"^(0^^',y) 



(7.24) 



l<fc<Af 



for 1 = 1,..., Af , where y e wp' and 

/ t)p)(x,y)d5x = 0. 
From this problem we can see 

Wo — 
/9n 



/ Mo^-^(x,y)dS'x = 
Jan Srix 



(7.25) 



(7.26) 



5f) 



0) 



^^'■'" "^^^^T^'^'^^l-es^^'^ "^^^ "*^°''^^^ 



for i = 1,... ,M. 



Before estimatingthe discrepancies in the boundary conditions, we note for 
y G wP , by Lemma ^ 



d 



Y: P('='°)(»7.)-^V,i?(^)(x,OW) 



l<fc<M 



9/1:, 



< const e, holds for x e dfle , 

(7.28) 



whereas the same lemma in conjunction with Lemma Wk), leads to 



'JV 



dn^ 



{ik,Vk) < const £2, forxe |J dJ^'^Udfl. (7.29) 



iyik 
l<i<M 
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Estimate of the right-hand side of ( 7.23 ) on 917. The definition oi J\f^^'^' in 
Section bl and the boundary condition (3.1) for R'-^\ give us 

Mo| ^{aA(^'''(^^, 7,^.) + (27rMo)-Mog(e- V - y|) + el'*''''(^,) • Vyi?("'(x, 0(^')}| 
d 



f^o 



9nx 



{-(2vr) ^(^j/ -^o^)log|x-y| 



■AjJ) 



~e{2TVfj.c 



~i-n(j,-f) 



V-^'^'i-qA- (x-O^^0|x-O 



»(j)|-2 



for X G 9f7, y G cj 



U) 



Using the asymptotics of h^' at infinity contained in Lemma 4p) , we obtain 



d 
9nx 



U) 



< const e , for X G dfl, y E u^ 



Then from Lemma [2] we have 

£Mo|^{2?(^''°)(l,)-V.i?(")(0(^),y) 

< const e^, X € dfl.y G ujy' . 



(7.30) 



l<fe<Af 



Therefore, this along with ( |7.28p , ( |7.29[ ) and ( |7.30| ) yield 



Mo 



9n, 



-(x,y) 



< const e , x e dfl, y G a; 



(i) 



(7.31) 



,(^) 



Estimate for the right-hand side of (7.24) on duje , i ^ 1, • • ■ , M . Consider 
first the situation when i — j. Since i?'-*'-' is smooth for x, y G 17, we can take 
the Taylor expansion of its first order derivatives about x = O'-'-' to obtain 



dRin) 



ap(j) 






yi^duJe 



'~a 



< const e, XG aa;p'',y Gwp'^ (7.32) 

where the boundary condition (3.9) for the dipole fields of the inclusion Scje 



was also used. Then, owing to Lemma [2] for i = j one has 



Oflx 



2^ V^^.^.;^(^fc)i 



l<fe<A/ 



dv^k.o) 



dn^ 



xe84^)--/^o^;^ 



a 



fc-'lx69uj(j) + 



Vxi?'"'(0(''\y 



< const e, X G 9w^-' -* , y G w^-' -* 
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This and (7.28), (7.29), (7.32) lead to the estimate 



M/.ir— (x,y)L ^. o)- -Moi^— (x,y) 



,(j)+ 



< const e, yew, 



0) 



(7.33) 
It remains to consider the case i ^ j. In this situation, we require Lemma [4J: 
to obtain 

9nx ^^''''" 



-($,,7,,) + (2^^o)-^^(log(e-V-y|)) 



d 






0<^'''(r7,)-{^Vx((2^Mo)"'log|x-0(^'r^)-^Vyii("'(x,0(^')}| 



< const e , X G Soj^'-' , y G w^^' i * 7^ J 



(7.34) 



The boundary conditions for the dipole fields (3.9) give 



(m/. -A'o) g^ (x,y) 

< Ia^/. - MO I |n^''^ • Vxi?(") (x, y) - V^i?(") (O^*) , y) | 

< const e, X e iSi^e , y € w^-'\ i ^ j ■ 



Then, Lemma [2] allows one to deduce 



aX'(i.o) 



aX>(j.o) 



A^/r 



(^.7) xPfli^w- -^^o^^- — «,-) 



E 

l<k<M 






artx ^^■'^l^eSo;^ 



Vxi?(")(0«,y) 

(7.35) 
•Vyi?(")(0(-'"),y) 



^^u■ 



dn^ 



aT)(k,0) 



< const e^, X e duj'^^\y e wp\ i 7^ j . 



Thus, from the preceding inequality and (7.28), (7.29), (7.34), (7.35), we have 
/^^. a^(^'y)lxeac.w--A*o^^(x,y)|^^g^(.,+ < const e, yew^J^^^J. 



We can conclude from this estimate, the conditions (7.25)-(7.27l, inequalities 
( |7.31| ), ( |7.33| , and Lemma [S] that 

M 

[)(^')(x,y)| <conste, x e |J c^^ U f]„y € c^p'^, j = 1, . . . , M . 
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f 7') 

Finally, the above estimate for 1)1 , j — l,...,Af and (7.22), combined with 
(7.121 complete the proof. D 



8 Numerical simulations 



In the current section, we implement the asymptotic formulae derived in Section 
[Tjfor Green's function N^ in numerical simulations. The numerical computations 
are carried out for the regular part i?^ of Green's function in {J^uji U V,^ for 
the transmission problem. In other words, let this regular part be given by the 
formula 



i?e(x,y) 



1 I xnAy) 



27r 



Mo 



M 

E 

(=1 



x^o (y) ^ 
^J■I, 



log|x-y| - Are(x,y) 



Then in accordance with the boundary value problem for iV^ given in Section 
[2] the function R^, which wc choose to consider for our numerical schemes, when 
y S il^ is a solution of the problem (7.3|-(7.8). 



8.1 Asymptotic formulae for R^ 

From formula ( |8.1[ ) and Theorem fl] we can immediately state the asymptotic 
formulae for the regular part R^ that will be used in the examples below. When 



y G ile we have from (7.1), that R^ admits the asymptotic representation 

M 



R, (x, y) = i?(") (x, y ) + 5] /j(^'°) {^^ , 77^. ) 



i=i 



M 



(i) 
When y e uje where j is fixed, j = I, . . . , M, R^ has the form 

i?,(x,y) = i?(^^)(x,y) + 4^^'(^^,,7,)+ E '^P''\^k,Vk) 



•2) 



k^3 
l<k<M 



+ {2n)-\^lo' - Ih!) loge - eV^=''\7^^) ■ Vyi?(")(x, O 



0)^ 



M 






i<fe<j\/ 



(8.3) 
Before proceeding with examples where these formulae will be implemented, 
we first discuss the numerical settings. 
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8.2 Numerical settings: Description of the geometry and 
physical parameters 

Let $7 be a disk of radius 150m, with centre at the origin, and occupied by a 
material with shear modulus fio = 5.6 x 10^"Nm~ , which is that of Cast Iron. 
We set the number of inclusions M ~ 6, and assume that the uj/ , j — I, . . . ,6, 
are circular. We summarize the data corresponding to the inclusions in Table 
1. 



Inclusion 


Centre 


Radius (m) 


Shear Modulus (xlOi"Nnr^) 


Material 


4^> 


(-90m, 40m) 


27 


2.6316 


Aluminum 


J'^ 


(-50m,-50m) 


24 


4.0741 


Copper 


J'^ 


(-30m, 10m) 


9 


7.7519 


Iron 


J/^ 


(20m, 70ni) 


19.5 


7.5188 


High Strength Alloy Steel 


J;'^ 


(50m, Cm) 


22.5 


8.0078 


Steel AISI 4348 


^(6) 


(70m, -80m) 


15 


9.0496 


Nimonic Alloy 90 



Table 1: Data for the inclusions 



,0) 



, J - 1, 



8.3 Model solutions used in the numerical simulations 
Neumann's function for the disk Vt 

In our examples, we need the Neumann function N^^i for a disk of radius R 
{R = 150m for our demonstrations), which is given by 



iV(")(x,y) 



log |x 



log 



i?2 



27r/io 27r/io 

and the regular part R^^' of this function is defined by the formula 

i?(^^) (x, y) = (2^/io)-i log |x - y| - N^'^^ (x, y) . 



The regular part of Green's function for the transmission 
problem in a plane with an inclusion at the origin 

Now, we state the form of the regular part hr^ of the Green's function for the 
transmission problem in the infinite plane with a circular inclusion at the origin 
of radius Qj . The solution is constructed using the involution procedure which 
is discussed in [5]. The representation of this function, is dependent on the 
position of the point ry -. When 77 ■ S Cw'^-'-' 



h^''''\i,,'n,) = - 




2TTflo /i/j 



l»7, 



for I, e Ca;(^) , 
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and 






For the case q^ e u}^^\ ^^ G Cuj'^^^ 
and for ?7j,$j G w*^^^ we have 



27r ^/^. + ^o V M/j Mo 



27r /i/ . + /io \ M/j 



loe 



IT^J 



Mo 



log flj 



The dipole fields for the circular inclusion in the infinite 
plane 

Here, we give the vector function V^^' whose components are the dipole fields 
for the circular inclusion of radius Cj in the infinite plane and 



2?(^''(l,) = Xcc.(.)(l,)2?(^'^°Hl,) +X.(.)(^,P^^''H^, 



We have 



pO-,o)(..) ^ i'l, f'oaU, ^^^ ^ ^_(,) ^ 



M/, + Mo II 



and 



v<^'''\i,) = 



M/, - f^o 
M/, + f^o 



^. , for 1^. e 



,0) 



8.4 Example 1 

The case of the force appUed outside the inclusions 

For our first example, we look at the case when y G il^. We therefore base our 



computations on the asymptotic formula (8.2) when comparing with those of 



COMSOL. The coordinates of the point force are given as y = (—10m, —80m). 
We plot the modulus of the gradient of the regular part R^ in Figure [2k) accord- 



ing to the analytical formulae (8.2 1. Figure 2b) is the same quantity computed 



using the method of finite elements in COMSOL. Both figures are very simi- 
lar, the maximum absolute error between these computations is 7.666 x 10~^^ 
occurring on the exterior boundary. 
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m^ 



i 



SftfT^ 






150 150 



a) 



150 150 



b) 



Figure 2: a) Computations based on asymptotic formula (8.2) and b) Numerical 
solution for the absolute value of the gradient of the regular part of Green's 
function for the transmission problem. Here y = (— lOni, —80m) and the mesh 
contains 44784 elements. The plots are practically indistinguishable. 



8.5 Example 2 

The case of the force positioned inside an inclusion 

In the second example, we aim to compare the computations produced by for- 
mula (8.3) with those generated by COMSOL. Now the point force is assumed 
to be situated at y = (60m, Om), in the inclusion Wg , containing the Steel AISI 
4340. Figure Isk), gives the surface plot for the modulus of the gradient of the 



regular part, provided by formula (8.3). The numerical solution given in COM- 
SOL is shown in Figure^). The maximum absolute error here is 7.98 x 10~^^, 

which occurs on the boundary of the inclusion uie ■ We conclude that the 
asymptotic formulae and numerical computations are in a good agreement with 
each other. 
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UodJlus of tte g^^die^tfort^ea^elytic^l x 




-150 -150 



a) 



[ 

■ -2.i 

■ -2 
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15D ' 150 



b) 



Figure 3: a) Computations based on asymptotic formula (8.3) and b) Numerical 
solution for the absolute value of the gradient of the regular part of Green's 
function for the transmission problem. Here y = (60m, Om) and lies in the Steel 
AISI 4340 inclusion. The mesh contains 44784 elements. There is once again a 
good similarity between the surface plots. 
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